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1 . 0  Research  Objectives 

1.1  A  New  Rapid  Acquisition  Technique  for  Direct 
Sequence  Spread  Spectrum  Communications 

The  rapid  acquisition  technique  described  here  can  be 
used  in  direct  sequence  spread  spectrum  systems.  The  tech¬ 
nique  employs  a  double  threshold  which  defines  when  a  decision 
can  be  made.  These  thresholds  change  at  each  examination 
instant.  Using  this  technique  a  significant  reduction  in 
the  acquisition  time  of  a  direct  sequence  spread  spectrum 
signal  is  obtained. 

1.2  A  New  Double  Threshold  Aquisition  Scheme  Applied 
to  the  Fading  Channel  in  Frequency  Hopping 
Spread  Spectrum 

This  technique  employs  the  double  threshold  procedure 
used  in  Sec.  1.1.  However,  the  design  and  results  obtained 
are  specific  for  frequency  hopping  in  a  fading  channel  such 


as  the  HF  channel. 


3) 


This  paper  describes  a  rapid  acquisition  technique  which  can  be  used  in 
direct  sequence  spread  spectrum  systems.  It  is  sequential  in  nature  and  is 
based  on  a  double  threshold  which  defines  when  a  decision  can  be  made.  The 
value  of  the  thresholds  change  at  every  examination  instant.  The  procedure 
makes  use  of  a  bound  to  the  partial  cross  correlation  function  of  PN  sequences 
which  insures  that  only  the  sequence  length  needs  to  be  given  as  far  as  the  PN 
sequence  particulars  are  concerned.  A  decision  is  made  only  if  the 
probability  of  error  associated  with  a  decision  at  that  point  is  equal  to 
or  better  than  a  prespecified  value.  Finally,  it  is  shown  that  most  of  the 
computational  difficulties  associated  with  the  bound  as  given  by  a  recursive 
algorithm  can  be  eliminated  when  using  a  low  order  polynominal  approximation 
to  the  bound.  However,  either  by  using  the  bound  or  its  polynominal 
approximation,  a  significant  improvement  in  the  acquisition  time  is  obtained. 

1.  INTRODUCTION 

The  problem  of  acquisition  of  a  spread  spectrum  signal  consists  of 
correctly  estimating  the  phase  of  a  received  maximal-length  shift  register 
sequence  (commonly  referred  to  as  a  PN  senuence)  of  arbitrary  length 
L  =  2N  -  1,  where  N  refers  to  the  number  of  delay  stages  in  the  generating 
shift  register.  The  acquisition  process  consists  of  repeatedly  comparing  the 
received  PH  sequence  /2P  g(t  -  j T  )  with  a  local  generated  PN  sequence  g(t 

j 

-  i Tc ) -  until  the  phases  iTc  and  jTc  are  found  to  be  equal.  At  that  point, 
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the  locally  generated  PM  sequence  is  said  to  be  synchronized  to  the  received 
PM  sequence  and  the  acquisition  process  is  ended.  Throughout  this  paper,  the 
carrier  phase  and  chip  timing  are  assumed  to  be  in  perfect  synchronism,  and 
the  only  interference  considered  is  additive  white  Gaussian  noise.  The 
in-phase/not-in-phase  decision  is  based  on  a  partial  correlation  process.  As 
is  well  known  [1],  the  partial  correlation  function  of  a  PN  sequence  is  not 
nearly  as  well  behaved  as  the  autocorrel ation  function,  it  being  a  function 
of  the  particular  starting  phases  i  and  j  ,  the  specific  feedback  connections 
of  the  shift  register,  and  the  length  of  the  correlation  yTc ,  where  y  equals 
the  number  of  chips  (an  integer)  and  Tc  equals  the  chip  duration.  To  circum¬ 
vent  all  these  difficulties,  various  bounds  have  been  found  ([2]-[4])  to 

the  partial  autocorrel ation  function,  and  one  of  these  bounds  is  used  in  the 
analysis  presented  in  the  next  section. 

In  particular,  a  sequential  estimation  technique  is  used  in  this  paper 
and  it  is  shown  in  Section  3  that  the  expected  time  to  make  a  decision  using 
the  sequential  procedure  is  upper  bounded  by  the  duration  of  time  needed  to 
make  a  decision  using  a  fixed  observation  interval  procedure  operating  with 
the  same  probability  of  error.  A  description  of  the  technique  as  well  as  the 
supporting  analysis  are  presented  in  the  next  section.  Numerical  results  are 
presented  in  Section  3,  and  Section  4  presents  the  conclusions. 

2.  ANALYSIS 

Reference  [3]  describes  a  method  for  generating  a  non-linear  sequence 
which  serves  as  an  upper  bound  to  the  partial  autocorrelation  function  of  a 
PN  sequence  (via  the  shi ft-and-add  property).  This  bound,  unfortunately,  is 
defined  via  an  algorithm  as  opposed  to  an  analytical  equation.  When  the 
incoming  sequence  is  correctly  aligned  with  the  reference  sequence,  the  output 
of  the  correlator  is  proportional  to  yT^  (see,  specifically,  (3)  below). 
Therefore,  for  simplicity  of  notation,  the  upper  bound  to  the  partial 
correlation  when  the  incoming  sequence  and  reference  sequence  are  not  in  phase 
will  be  written  as  yTc(l  -  G(y)).  This  then  results  in  yG(y)  as  being  the 
(normalized)  difference  between  the  in-phase  and  the  not-in-phase  correlation 
curves. 
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In  addition  to  presenting  precise  results  based  upon  the  upper  bound, 
approximate  results  will  also  be  presented.  These  latt'r  results  arise  from 
curve  fitting  the  exact  bound  with  a  low  order  polynomial  in  y.  Specifically, 
it  can  be  shown  that 

yC  1  -  G(y)]  =  yU  -  pi-  (1) 

where  L  is  the  period  (i.e.,  the  length)  of  the  Pil  sequence. 

In  Section  3,  numerical  results  using  the  approximate  bound  will  be  compared 
to  those  obtained  using  the  exact  bound. 

Referring  now  to  Figure  1,  the  procedure  works  in  the  following  manner: 

An  integration  time  y^  (in  integer  number  of  chips)  is  established  so  that  a 

given  probability  of  incorrect  decision  Pg  can  be  achieved  with  that 

integration  time.  Assuming  that  Yj  >  N,  where  N  is  the  number  of  stages  of 

the  shift  register  which  generates  the  spreading  sequence,  then  if  any 

integration  time  y2  such  that  N  <  y2  <  y^  is  chosen,  Pg  cannot  necessarily  be 

achieved.  However,  by  allowing  an  uncertainty  region  such  that  if  the  test 

statistic  falls  in  that  uncertainty  region  then  no  decision  is  made,  but  rather 

the  integration  time  is  increased  by,  say,  one  chio,  the  test  can  be  reapplied  with  the 

increased  integration  time  to  again  see  if  Pg  can  be  achieved.  This  process 

is  repeated  continually  until  an  integration  time  of  y^  is  used,  at  which 

point  Pfi  is  guaranteed  to  be  achieved.  This  will  become  more  clear  as  the 

details  of  the  procedure  are  described  below. 

The  output  of  the  correlator  in  Figure  1  is  qiven  by 

/yj 

o  C  Y&l  S(t  -  JTC>  *  nwU>]  sit  -  >TC)  dt,  (2) 

which  consists  of  a  signal  term  sQ(yTc)  and  a  noise  term  nQ( yTc ) .  The  noise 
term  can  easily  be  shown  to  be  Gaussian  with  zero  mean  and  variance  02(yT  )  = 

7yTc’  where  \  is  the  two-sided  power  spectral  density  of  the  input  noise. CThe 
signal  term  sq(yTc)  is  given  by 
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r^c 

So(yTc)  J  /2Fl  9(t  '  jTc}  *  9<t  "  iTc)dl 


1=  /2fVTc 


j  f  i  =  j 


S  ^2PsYTcU-GlYn  if  i  i  j 
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where  the  lower  term  on  the  rhs  arises  from  the  bound  in  (1).  Note  that  all 
the  following  equations  which  depend  upon  this  bound  are  themselves  bounds. 
For  simplicity,  however,  they  are  wr.itten  with  equality  signs. 

The  synchronization/no-synchronization  decision  is  seen  only  to  involve 
the  choice  of  which  signal  term.’s  (YT  ),  is  present  at  the  correlator's 
output.  Due  to  the  presence  of  noise,  this  decision  can  only  be  made  with  a 
specified  accuracy.  If,  as  shown  in  Figure  2,  we  wait  for  a  time  T  *  yjTc 

before  making  a  decision,  then  we  can  define  a  threshold  such  that  if  the 

correlator's  output  falls  above  the  threshold,  we  decide  synchronization  is 

attained,  and  if  it  falls  below  the  threshold,  we  decide  the  two  PN  sequences 

are  not  synchronized.  The  probability  of  an  erroneous  decision,  given  we  are 
not  looking  at  the  correct  phase  position,  is  equal  to 


Pe  =  P[Wc>  >  V  =  P["o(Vc)  >  (4) 

where  ej  is  shown  in  Figure  2  and  VT,  the  threshold,  is  given  below  [see 

(7)3* 

Note  that  both  in  Figure  2  and  from  Equation  (4),  for  simplicity,  the 
decision  boundary  has  been  located  halfway  between  the  two  signals.  There¬ 
fore, 
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The  probability  of  error  of  Equation  (4)  then  becomes 

.%  ■  "K'nV  »  7  Vc  G(>in-7erfc 

where  E$  -  y.Ec,  £c  is  the  energy  per  chip  and  equals  PsTc,  and 

oo  2 

erfc(x)  -  —  f  e  y  dy. 

J  x 

Finally,  Vj(tiTc)  is  given  by 

VnV  ■  nTc  -  *i  *  w;  nt [1  -  r  Gln)] 


(6) 


(7) 


The  following  discussion  will  show  how  the  acquisition  process  can  be 
improved  by  shortening  the  acquisition  time. 

The  Decision  Process  prior  to  y.  =  y^:  At  times  it  may  not  be  necessary 
to  wait  the  full  chips  to  make  a  decision.  The  only  requirement  is  that 
at  any  time  a  decision  is  made,  the  probability  of  error,  Pg,  be  kept  con¬ 
stant.  If  any  attempt  at  raking  a  decision  after,  say,  y2  chips  is  made, 
then  the  threshold  voltage  must  be  set  such  that 

Pe  =  P[no(Y2Tc)  >  c23  =7  erfc 


^o^V 


(8) 


Note  that  Equation  (8)  is  identical  to  Equation  (4),  with  the  exception  of 
the  pertinent  parameters  (i.e.,  the  decision  mechanism  is  identical).  But  if 
the  probability  of  error,  P  ,  is  to  be  kept  constant,  then 


at  t 


y?' 


(9) 
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This  implies  that  the  two  arguments  of  the  erfc  function  must  be  equal.  Hence 


or 


^0(Vc>  ^Oo(  y2Tc} 


(10) 
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e2 =  ci  =  (nTc)  GlV 


J1 

^Yi 


Equation  (10)  defines  the  parameter  e£  such  that  at  any  time  T  =  Y2TC  a 
decision  can  be  attempted.  Figure  2  shows  qualitatively  the  decision  regions 
at  any  given  ^  <  Y1 ' 

If  at  time  T  =  y2tc  the  correlator's  output  voltage  equals  VQ,  the 
following  decisions  can  be  made  with  a  given  probability  of  error  Pg: 

1)  Synchronization  has  been  attained  (i  =  j)  if 

*o  >  VS  ■  y2Tc(1  -  G(-'Z)>  *  c2 

2)  Synchronization  has  not  been  attained  {i  f  j)  if 


Vo  <  Vn  =  /2FV2Tc 


c2 


3)  No  decision  can  be  made;  continue  the  correlation  if 


Thus,  Vn( y2>  and  V$(y2)  become 


W  =  ^  '2Tc  -  e2 
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and 


V(y2)  =  22T  Tc  ,2(l  -  S(t2>)  =VWcC2vf  (l  -  «r2>)  *  Gty^)  . 


(12) 


respectively. 


From  Figure  2,  the  probability  of  jto  decision  can  be  seen  to  equal 


P  (no  decision)  -  Pnd  =  P[Vn  <  sQ  +  nQ  <  vs^ 


(13) 


With  Vp  and  V$  given  in  Equations  (11)  and  (12),  respectively,  the  probabil¬ 
ity  of  no  decision  at  time  Y-Tc  becomes 


Of  more  interest  than  (14)  is  the  conditional  probability  of  no  decision 
at  the  jth  step,  given  that  no  decision  had  been  made  at  the  j-lst  step.  For 
ease  of  notation,  let  i  -  j-1,  and  denote  by  fg  and  fp  the  value  of  the  in¬ 
sync  and  not-in-sync  curve  of  Figure  2,  respectively,  Jat  the  point  Yj-  De_ 
noting  the  probability  referred  to  above  as  P-j,  we  have 
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where 


V  . 
SJ 


VV 


un 


and 


y*T 

rij  =  J  j  Cnw( t )  g  ( t  -  kTc)dt,  (16) 

where  k  corresponds  to  the  specific  phase  position  of  the  spreading  sequence 
that  is  currently  being  examined. 

The  noise  output  at  time  t,  is. given  by 

«l 

"j  ■  ni  * 4n 

where  An  is  the  noise  contribution  in  the  interval  At  (e.o,,  At  =  Tc) , 
and  the  threshold  voltages  are  given  by 


and 


Vsj  "  fnj  +  Ej 


Therefore,  (15)  can  be  rewritten  as 


P,j  •  P  C-Cj  -«,<»»<  (f„j  -  fsj)  *  ‘j  -  "j  |  <  Ui  *  ni  <  vsi] 

■  P  {-'j  -  "i  <in  '  (frj  -  fsj*  *  ej  -  V  -  £i  <  "i  <  (f„i  -  fsi>  *  pn 
’H l-'i <  "i  <  rf„7--fsi-r.-7-|  v-t 


If,  for  simplicity,  0,  S  Gt^.),  then  the  numerator  of  (19)  can  be  shown 

J  J 


to  reduce  to 
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(19) 
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be  the  event  {  making  a  decision  for  the  first  time  at 
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(22) 


where  is  given  by  (19)  and  Pnd  (yj  is  6iVerlty  (H,).  Illicit  in 
(22)  is  that  the  first  time  a  decision  is  attempted  is  y.  =  N. 
Finally,  it  is  of  interest  to  determine  the  average  time  needed  to  make 


a  decision.  To  evaluate  this  quantity,  the  Markov  chain  model  described  in 
[5]  can  be  used.  The  transition  diagram  is  shown  in  Figure  3,  where  states 
N+l,  N+2,  ....  (yj-1)  indicate  both  a  possible  path  to  state  D,  the  decision 
state  (absorbing  state),  and  also  a  possible  path  to  the  next  higher  state. 
On  the  other  hand,  when  state  y1  is  entered,  the  next  state  has  to  be  D. 


The  canonical  form  of  the  transition  matrix  (see  [5])  for  this  scheme  is 
given  by 
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where  the  partitions  5,  0,  R,  and  Q  are  defined  in  [’I.  If  we  define  the 
matrix  N  -  [I  -  Q]  where  I  is  an  identity  matrix,  it  is  shown  in  [5]  that 

r  =  n  1  (24) 

L  L. 

where  ^  is  a  column  vector  whose  jc  element  is  the  mean  time  to  a  decision 

given  the  process  started  in  state  j,  and  \  is  a  column  vector  whose  elements 
are  the  duration  of  time  needed  to  make  each  transition. 


For  this  particular  acquisition  scheme, 


C  = 


(25) 


Also,  it  is  only  the  first  entry  of  the  vector  t  that  is  of  interest, 
since  we  are  interested  in  the  average  time  to  a  decision  given  that  we  start 
in  state  [>(.  From  {  23) .  (24  ),  (25  ),  and  the  definition  of  F),  we  can  show  that 
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T1  =  *N+  PN,N  +  1  +  PH  ,N  +  1PN  +  1  .N-t-2  +  ...+  P;.'  , i!  + 1  ?N  + 1  , li • 
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Notice  that  is  upper  bounded  by  YjTc>  which  would  be  the  value  of  the 
integration  time  if  a  fixed  interval  was  used  to  make  the  decision. 

3.  NUMERICAL  RESULTS 

In  this  section,  the  key  equations  describing  the  performance  of  the 
system  will  be  evaluated  numerically  to  illustrate  how  the  acquisition  scheme 
behaves.  The  parameters  are  chosen  arbitrarily  to  be  indicative  of  realistic 
values  on  one  hand,  but  to  not  result  in  excessive  computation  time  on  the 
other  hand.  In  particular,  there  are  six  figures  presented,  each  one 
containing  two  curves  graphed  on  a  common  set  of  axes.  These  curves  represent 
the  probability  of  making  a  decision  within  y  chips.  Results  are  presented 
for  the  exact  bound  as  well  as  for  the  approximation  to  the  exact  bound.  The 
curves  plot  the  probability  of  a  decision  along  the  ordinate;  the  number  of 
chips  within  which  the  decision  is  made  >s  the  abscissa.  For  each  curve  we 
have  set  the  length  of  the  PN  register  to  be  10  so  that  the  period  of  the 
sequence  is  L  =  1023.  The  curves  are  then  plotted  for  a  chip  energy-to-noise 
density  ratio  of  0  dB  to  -20  d!3  and  a  probability  of  error  of  either  0.1  or 
0.01.  However,  using  Eq.(22),any  other  set  of  parameters  may  be  employed. 

It  can  be  seen,  for  example  from  Figure  6,  that  if  the  ratio  energy- 
per-chip-to-noise  density  E^/q  is  -20  dB  and  the  probability  of  error  is  .1, 
then  for  a  probability  of  decision  of  0.99,  the  two-threshold  sequential 
procedure  will  result  in  detection  within  140  chips,  compared  to  a  single 
threshold  fixed  dwell  time  scheme  which  requires  an  interval  of  560  chips 
(with,  of  course,  a  probabilty  of  decision  of  1.0).  If  the  probability  of 
error  is  decreased  to  0.01,  it  is  seen  from  Figure  9  that  the  two-level 
technique  will  result  in  detection  within  520  chips  whereas  the  single 
threshold  fixed  dwell  scheme  requires  820  chips. 


If  the  approximation  curves  shown  in  Figures  4  through  9  are  compared  to 
the  corresponding  exact  curves,  it  will  be  seen  that  at  low  values  of  tjn  and/ 
or  higher  values  of  Pe,  the  approximation  yields  fairly  accurate  results.  On 
the  other  hand,  for  lower  values  of  Pe  and/or  higher  values  of  E  / n •  the 
approximation  is  rather  loose.  However,  the  simplicity  of  evaluating  the 
approximation  together  with  the  incite  one  can  obtain  by  having  a  closed-form 
expression  (as  indeed  is  the  case  with  the  approximation)  make  it  a  useful 
tool  in  this  type  of  design. 

Finally,  it  is  seen  from  some  of  the  figures  (e.git  Figure  4)  that  the 
upper  curve,  which  corresponds  to- the  exact  upper  bound,  has  a  very  irregular 
shape.  This  is  due  to  jumps  in  the  upper  bound  itself  and  to  the  discrete 
nature  of  the  computation  (that  is,  due  to  the  fact  that  each  additional 
increment  of  integration  time  is  of  length  one  chip). 

4.  CONCLUSION 

In  this  paper,  we  have  analyzed  a  rapid  acquisition  technique  for  DS 
spread  spectrum  communications.  The  technique  employed  was  sequential 
estimation,  and  used  two  thresholds  at  each  step  in  the  procedure  to  determine 
whether  or  not  synchronization  was  to  be  declared;  the  probability  of  error  at 
each  step  in  the  process  was  always  held  constant. 

Numerical  results  were  presented  which  illustrated  the  performance  of 
this  approach.  We  have  determined  the  probability  of  making  a  decision  within 
y  chips  and  have  proven  that  the  average  time  to  making  a  decision  using  ‘this 
technique  is  always  upper  bounded  by  the  time  needed  to  make  a  decision  using 
a  fixed  interval  acquisition  procedure. 

Specifically,  Figures  4-9  illustrated  the  actual  savings  in  acquisition 
time  which  may  be  expected  for  a  sequence  of  length  1023,  a  chip-to-noise 
density  ratio  of  0  dB  to  -20  dB,  and  an  acceptable  Pe  of  either  0.1  or  0.01. 
These  results  showed  that  savings  of  40%  to  75%  in  the  acquisition  time  are 
entirely  feasible  if  the  double  threshold  scheme  is  used  instead  of  a  single 
threshold  procedure. 
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ABSTRACT 


In  this  paper  a  double  threshold  acquisition  scheme 
is  applied  to  the  fading  channel  in  frequency  hopping  spread 
spectrum  in  order  to  minimize  the  mean  acquisition  time.  The 
model  of  the  fading  channel  is  presented.  The  acquisition 
strategy  is  described.  The  equations  for  the  two  thresholds 
are  derived.  A  computer  simulation  program  is  used  to 
determine  the  mean  and  the  variance  of  the  acquisition  time. 
The  results  are  analyzed  and  compared  to  the  nonfading 
channel  and  also  compared  to  the  single  threshold  scheme  in 
the  fading  channel. 


DOUBLE  THRESHOLD  ACQUISITION  SCHEME  APPLIED  TO  THE 


FADING  CHANNEL  IN  FREQUENCY  HOPPING  SPREAD  SPECTRUM 


1.  INTRODUCTION 

In  FH  communication  system  the  transmitting  carrier 
frequency  is  controlled  by  a  pseudonoise  sequence  generated 
by  linear  feedback  shift  registers.  The  state  of  the  shift 
registers  controls  a  frequency  synthesizer  which  provides 
the  transmitting  carrier  frequencies.  In  the  receiver  there 
must  be  such  a  PN  sequence  replica  which  controls  the  local 
carrier.  The  two  PN  sequence  should  be  synchronized  such 
that  the  local  carrier  frequency  synthesizer  provides 
exactly  the  same  hopping  frequency  and  hops  at  the  same  rate 
as  the  transmitter  in  order  to  dehop  the  signal  and  recover 
the  data  from  demodulation  such  as  FSK .  Before  the  system  is 
synchronized  no  data  can  be  sent.  Therefore,  the  mean 
acquisition  time  is  greatly  concerned.  The  purpose  of  this 
paper  is  to  minimize  the  mean  acquisition  time.  Based  on  the 
work  that  we  have  done  to  the  nonfading  channel  r  reference 
11  a  double  threshold  acquisition  scheme  is  also  applied  to 
the  fading  channel.  A  computer  simulation  program  is  used  to 
determine  the  mean  acquisition  time  and  the  variance  with 
respect  to  the  other  parameters.  Tn  this  paper  the 
acquisition  strategy  is  described.  The  faded  signal  has  been 
studied.  The  probability  density  function  of  the  signal 


envelope  square  is  given.  The  equation  of  the  two  thresholds 
are  derived.  The  computer  simulation  is  done  statistically 
for  each  of  the  specific  parameter.  The  results  are  analyzed 
and  compared  to  the  nonfading  channel  and  also  compared  to 
the  single  threshold  scheme  in  the  fading  channel. 


2.  ACQUISITION  STRATEGY 

The  tracking  loop  works  only  when  the  alignment 
between  the  transmitter  and  the  receiver  is  greater  than 
half  of  the  hop  duration;  therefore,  the  acquisition 
technique  includes  two  stages,  pre-acq  and  true-  acq,  to 
gaurantee  this  requirement. 

Acquisition  strategy  of  FH  signal  is  the  following: 
Refer  to  Fig.  1  and  Fig.  2. 

The  incoming  signal  having  frequency  and  phase  0  is 

embedded  in  noise.  The  local  frequency  synthesizer  is 

multiplied  by  the  incoming  signal  and  the  product  passes 

through  a  BPF  which  eliminates  the  sum  frequency  components 

and  presents  the  difference  component  and  the  noise 

1  1 

term  to  the  integrator.  and  are  selected  such  that  the 

VO  -j  Uu  _ 

x  J 

output  of  the  integrator 
I -  ,  rT 

J^s  J  COS  (ub  )  t  dt  —  0  ,  tor  bl  . 

0  x  J  i  j  (P-1) 


ie ,  (f . -f  . ) T 

i  D  C 


and 


integer 


COS  (<jJ.  -eg,  it  dt  =  rT  p 
i  j  c  ,  s 


(2-2) 
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The  noise  has  a  mean  value  of  0  and  a  variance  of  n /?. 

After  integration  the  noise  term  present  at  the  output  of 

2 

the  integrator  has  a  mean  value  of  0  and  the  variance  0  is 
proportional  to  hrT^/2. 

The  signal  has  two  components,  in  phase  and 
quadrature.  We  then  take  the  envelope  (or  envelope  square). 
The  envelope  square  is  compared  with  two  thresholds  THA 
(lower  threshold)  and  THB  (upper  threshold).  Refer  to  Fig  3, 
at  time  t^  the  output  of  the  integrator  is  less  than  THA, 
we  say  that  the  incoming  signal  is  not  synchronized  with  the 
local  frequency  synthesize.  The  comparator  outputs  a  signal 
which  controls  the  local  PN  generator  to  shift  to  the  next 
state,  thus  the  frequency  synthesizer  also  hops  to  the  next 
frequency  f^.  at  time  t?  is  greater  than  THB,  we  say  that 
pre-acquisition  has  occurred.  The  frequency  synthesizer  and 
PN  generator  remain  in  the  old  state  until  the  incominq 
signal  hops  to  the  next  frequency  f?.  At  time  t^  v0  aqain 
drops  below  THA,  the  comparator  declares  that  true 
acquisition  has  occurred,  the  PN  generator  shifts  to  next 
state  and  the  frequency  synthesizer  hops  to  next  frequency 
f 2 •  In  this  case  the  receiver  is  at  most  one  chip  (rt^) 
later  than  the  transmitter;  that  is,  the  alignment  between 
the  incoming  signal  and  the  local  frequency  synthesizer  is 
longer  than  half  of  the  hop  duration.  At  this  moment  the 
tracking  loop  starts  to  work  which  brinqs  the  local  hopper 
completely  in  phase  with  the  incominq  signal  and  thus,  fine 
acquisition  is  fulfilled. 
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If  is  between  THA  end  THR  ,  no  decision  is  mede  and 
more  samples  must  be  taken.  At  the  end  of  the  second  r?c 
another  comparision  is  made.  ...etc.  This  procedure  is 
repeated  until  a  decision  is  made. 

In  our  simulation  rT^=  T^  fone  bit  duration). 

3.  PROBABILITY  DENSITY  FL'T^TTON  OF  THE  SIGNAL  ENVELOPE 
IN  THE  FADING  CHANNEL  AND  THE  TWO  THRESHOLDS 


The  mean  value  of  the  coefficients  of  the  fading 
channel  is  zero  so  we  say  that  the  fading  channel  without 

ie 

specular  component  is  Rayleigh  channel. 

In  general,  if  X  and  Y  are  Gaussian  random  variable 

2 

with  variance  0  . 


(3-1  ) 


2  2 

f  (Z)  =  e_z  / 2 6  U(z)  - Rayleigh  Density 

<5  2  (3-2) 

U ( z )  =  unit  step  function 

For  convenience  we  consider  the  envelope  sauared.  We  define 

g 

W  ~  7  ,  W  =  envelope  square 

The  probability  density  function  of  W  is; 


f  (W) 


1  -W/2  62 

262 


U(W) 


>3-3) 


*:  See  Appendix  A 


The  probability  tna t  W  is  between  a  anr.  is: 


P(c/^W  <  (3  )  =  ’  f(W)dW  =  e'U/262  -  e-/5/2!5' 


n--2) 


Before  synchronisation  we  have  noise  only. 


variance  equals 


When  the  system  is  going  to 


synchronize,  we  have  signal  and  noise.  The  variance  equals 


2  *  * 


Refer  to  Fig. when  the  system  is  not  synchronized, 


for  each  comparision  we  make  an  error  if  the  noise  is 
greater  than  TH3. 


PFA=  p(noise  £  THB)  =  f  (W)  cW 

TiiB 


n-5) 


When  the  system  is  synchronized,  we  do  not  acknowledge 
the  synchronization  if  the  signal  +  noise  is  less  than  THA. 


PD  “  ?  (signal  +  noise  4  THA)  = 


f  (Vi)  aW 


=  1  -  e 


-THA/  2  5  2 

SN 


(7-5) 


From  eqs  (2-5)  and  (2-5)  we  have 


TKA  *  -2(5  2 


In  (  1  -  p_ 


THB  =  -2(5  .  in 

«  iA 


(2-7) 


.2-8) 


Fig.  4  shows  the  probability  density  function  of 
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envelope  square,  P^  and  PD<  Fiq.5  and  Fig.B  give  THA  and 
THB  versus  number  of  bits  being  integrated  for  the  different 
signal  to  noise  ratio  and  the  different  P^. 

In  order  to  minimize  the  probability  of  false 
acquisition  P0.  will  be  less  than  or  equals  to  0.0001.  P_ 

r  A  D 

affects  the  mean  acquisition  time  significantly.  This  can 
happen  in  two  ways: 

Refer  to  Fig  a, 

a)  With  a  small  PD  most  of  the  noise  will  be  above  THA  and 
we  always  cannot  make  a  decision  which  increases  the  mean 
acquisition  time  . 

b)  If  PD  is  large  the  signal  +  noise  will  also  lie  below  THA 
and  it  will  be  falsely  dismissed. 

Both  of  the  above  situation  will  increase  the  mean 
acquisition  time. 

4.  COMPUTER  SIMULATION  RESULTS  -  MEAN  ACQUISITION 

TIME  AND  VARIANCE  VERSUS  Eb/n  FOR  DIFFERENT  PR. 

The  normalizd  mean  acquisition  time  versus  E^/  n  for 
the  different  PR  are  shown  in  Fig. 7  and  Fig.fi.  Fiq.7  is  for 
30  hops/S  and  Fig.B  is  for  100  hops/S.  From  these  figures  we 
observe  that  the  mean  acquisition  time  doesn't  depend  on  the 
hop  rate  and  mostly  depends  on  the  energy. 

Comparing  with  the  nonfadinq  channel  f  reference  11  we 
observe  that  when  Eb/n=3db,  fading  increases  the  mean 
acquisition  time  ??%  only  because  comparing  to  the  noise 
fading  is  not  important.  When  Eb/n=70db,  fading  doesn't 


affect  the  mean  acquisition  time  at  all.  The  mean 
acquisition  time  is  close  to  the  minimum  :  ■  value 

it  it  1c 

which  equals 

L=  code  length,  in  our  simulation  L=hll 
(L/2)  *1 .  1  =  f  S 3  1  / 7 )  *  1  .  1  =  ?. a ]  bits 

When  E^/n  equals  iSdb,  equals  0.7  and  0 .  S , (points  A 
and  B  in  Fig. 7)  and  also  refer  to  Fig.?;  we  observe  that  THA 
and  THB  intersect  earlier  than  or  at  the  vicinity  of  one 
bit.  Considering  the  acquisition  strategy,  refering  to  Fig.? 
,  we  would  find  that  acquisition  only  depends  on  THB.  Tn 
that  case  the  double  threshold  scheme  becomes  single 
threshold  scheme  and  the  mean  acquisition  time  goes  high. 
Point  C  ( PD=0 . 1 )  is  still  the  case  of  double  threshold 
scheme.  Hence,  this  is  the  improvement  of  the  double 
threshold  scheme  comparing  to  the  single  threshold  scheme. 
Again  at  points  D  and  E  double  threshold  scheme  becomes 
single  threshold  scheme  and  we  may  reduce  the  PD  to,  say, 
0.02  to  find  another  double  threshold  scheme.  Point  F  is  the 
improvement . 

The  variances  of  the  acquisition  time  are  listed  in 
Table  1  and  Table  2. 


***:  See  Appendix  C 
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APPENDIX  A 


MODEL  OF  THE  FADING  CHANNEL 


The  fading  channel  is  represented  by  two  tapped  delay 
lines,  one  representing  the  in  phase  component  and  another 
representing  the  quadrature  component.  To  represent  a  HF 
channel  the  taps  are  spaced  about  1ms  apart.  There  are 
either  1  or  4  taps  employed.  The  signal  samples  from  those 
taps  are  multiplied  by  independent  Gaussian  noise  processed 
as  shown  in  Fig. 9.  The  output  signal  from  the  fading  channel 
is  then  represented  by  two  equations. 


1  NSC  DSP(j)*CP(j)-OSQ(j)*CQ(j) 
j=l,TDS,— 

(A-l) 

£  __  DSQ  (j  )  *CP  (j  )+DSP  (j  )  ,VCQ  (j  ) 

j=l,TDS,-7p- 

(A-2) 

In  equation  (A-l)  and  (A-2) 

FP:  Fading  in  phase  sample 
FQ:  Fading  quadrature  sample 

TDS :  Total  number  of  delayed  sample,  in  phase  or  quadrature 

i 

TDS  =  (  (NEAPS- 1)  *  NSC  )  +  1  (A-3) 

NTAPS :  Total  number  of  taps,  in  phase  or  quadrature 

NTAPS  =  CHIP  RATE  ‘CHANNEL  SPACING  *  SPREAD  TIME  *7  fA-4) 

NSC:  Number  of  samples/chip 

DSP(J):  Individual  in  phase  delayed  sample 

DSOM):  Individual  quadrature  delayed  sample 


v 

U  , 
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CPM):  Individual  in  phase  coefficient 
COM):  Individual  quadrature  coefficient 
RP:  Transmitting  in  phase  signal  sample 
RQ:  Transmitting  quadrature  signal  sample 

CP  and  CO  are  random  variables  which  represent  the 
character i st i c  of  the  fading  channel.  CP  and  CO  equals  a 
Gaussian  random  number  multiplied  by  its  standard  deviation 
and  then  passes  through  a  low  pass  filter  whose  cutoff 
frequency  equals  the  fading  rate.  The  data  passing  through 
the  low  pass  filter  will  lose  energy  and  therefore  its 
variance  is  adjusted  accordingly.  The  adjustment  is 
proportional  to  the  sampling  frequency  and  inversely 
proportional  to  the  fading  rate.  The  constant  of 
proportionality  is  about  1/M.  The  variance  of  CP  and  CQ 
equals  the  total  normalized  power  divided  by  the  number  of 
taps  and  multiplied  by  the  factor  of  the  adjustment. 


Variance  before  filtering  -  *  5SgS-H!£  *  * 


(A-S) 


Variance  after  filtering  =  g^APS 


(A-fi) 

In  the  case  of  FSK,  as  shown  in  Fiq.lh,  R0=0  ,  Fqs 
(A-l)  and  (A-7)  become  Fqs  (A-7)  and  (A-R). 


FP  = 


j=l ,TDS , 


NSC 


DSP (j )  *  CP(j) 


(A-7) 


FQ  =  L  Nsc 


DSP(j)  *  CO(j) 


(A-8) 
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APPENDTX  B 

SIGNAL  IN  THE  FADING  CHANNEL  AND  THE  VARIANCE  OF  THE 
RECEIVED  SIGNAL  AT  THE  OUTPUT  OP  THE  INTEGRATOR 


The  fading  channel  is  represented  by  two  delay  lines. 
When  the  transmitter  hops  to  a  new  frequency,  the  signal  is 
passed  through  a  new  fading  channel.  Consider  the  in  phase 
component  only.  (The  quadrature  component  is  the 
same) .Refering  to  Fig. in,  for  a  2-tap  model  we  see  that  the 
first  to  the  fifth  fading  samples  are: 


FP(1)  =  JI  *  CP(1) 


FP(5)  =  JI  *  0P(1) 


=  JI  *  _J_ 

J2NTAPS 


=  0.707 


(B-l) 


The  sixth  to  tenth  fading  samples  are 


FP(6)  =  Z  JI  CP(j) 
3=1 


(B-2) 


Because  CP's  are  Gaussian  random  variables  with  a 
variance  of  1/2NTAPS  and  a  frequency  bandwidth  of  H^.,  the 
fading  rate,  the  summation  among  the  different  taps  is 
calculated  according  to  the  power.  Hence, 


FP(6)  =  JI  Z.  CP(j) 
j=l 


=  JI  *  JI  * 


'2NTAPS 


(B-3) 
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For  FP(ll)  to  FP{100),  since  ail  the  delay  iines  are 
filled  with  signal,  the  magnitude  of  FP  will  be  3. 

At  the  time  of  the  101th  sample,  the  transmitter  hops 
to  the  next  frequency  and  the  delay  lines  are  gradually 
filled  by  new  values.  Nevertheless,  some  delay  lines  are 
still  filled  with  the  old  frequency  signal,  however,  the 
magnitude  of  the  old  frequency  signal  decreases  gradually 
according  to  the  same  rule. 

For  the  )  01th  to  105th  sasmples 

FP(lOl)  -  J?  CP(2)  =  p  *  Jg-L  -  0.707  (B-4) 


For  a  d-tap  model,  NTAPS=d  ,  similarly, 

FP (1 )  to  FP (5) =0 . 5 
FP  ( 0 )  to  FP  (]0)  =0.70-7 
FP (11)  to  FP(15)=0.055 
FP (15)  to  FP ( 100) =] 

FP (101)  to  FP (105) =0.005 
FP (105)  to  FP (1)0) =0 .707 
FP(lll)  to  FP ( 1 1 5 ) =0 . 5 

The  signal  variation  with  time  in  the  fading  channel  is 
shown  in  Fig. 11  to  Fig. 33  for  the  different  hopping  rates. 

When  the  system  is  not  synchronized,  the  received 
detected  output  consists  of  noise  only.  When  the  receiver 
hops  to  the  transmitting  frequency,  there  are  ten  possible 
situations  if  we  have  ten  bits  in  one  hop.  The  receiver  can 
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only  hop  at  the  end  of  the  bit  time.  If  the  transmitter  in. 
the  receiver  hop  simultaneously,  there  is  no  signal  in  the 
channel;  otherwise,  some  samples  are  already  in  the  fading 

t 

channel.  Now  we  consider  the  situation  that  both  the 
transmitter  and  the  receiver  hop  simultaneously.  At  the  end 
of  the  first  bit  the  receiver  receives  ten  samples.  The 
first  five  samples  are  correlated  because  they  have  the  same 
CP(1).  Due  to  the  low  fading  rate  these  five  samples  add 
almost  linearly.  (The  precision  is  up  to  the  third  digits). 
Sample  fi  to  sample  10  also  add  linearly.  But  sample  1  to 
sample  5  and  sample  6  to  sample  10  only  add  partly  linearly 
because  they  have  the  common  coefficient  CP(1)  but  samples  ] 
to  5  don't  have  CP(2).  For  simplicity  we  suppose  all  the 
samples  add  linearly.  Of  course,  this  is  an  approximation. 
But  the  influence  of  such  approximation  will  be  partly 
cancelled  if  the  signal  is  already  in  the  fading  channel 
when  the  receiver  hops  to  the  transmitting  frequency. 
Therefore,  for  a  two-tap  model  at  the  end  of  the  first  bit 
the  standard  deviation  of  the  received  signal  at  the  output 
of  the  integrator  is: 

^sl  “  *  ft s  ] 2NTAPS  *  ~=-  +  [P~  .  *  MSC  j  p  =2 

2  J  s  J 2NTAPS  2  s 
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From  bit  2  to  bit  10,  becauses  all  the  delay  lines  are 
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lied  with  signal;  hence, 
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for  bit  11,  the  signal  in  the  delay  lines  gradually  move 
out;  hence. 
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(3-7) 


For  a  four-tap  model: 
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6S3  =  6s2  +  NSC 
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As  for  the  noise: 

u  L. 

At  the  end  of  the  Mc  bit,  the  standard  deviation  of  the 
noise  is: 


WNL 

^NM  =  VM  *  NSC  *NCB  *  10  20 
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WNL:  Eb  /  in  db 

NCB:  number  of  chips/bit,  which  equals  1. 

At  the  end  of  the  Mc  bit,  the  standard  deviation  of  the 
signal  plus  noise  is: 
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APPENDIX  C 

THE  MINIMUM  POSSIBLE  ACQUISITION  TIME 


If  there  is  no  noise  the  minimum  possible  acquisition 
time  equals  (L/?)*1.1 
L  =  code  length 

The  average  difference  between  the  transmitter  and  the 
receiver  is  half  of  the  code  length  L/?.  The  factor  1.1  is 
due  to  the  hopping  rate  of  the  transmitter  during 
acquisition.  If  there  is  no  noise  so  that  the  receiver  can 
make  a  decision  at  the  end  of  each  bit,  it  needs  L/2  bit 
time  to  catch  up  the  transmitter.  But  in  this  time  interval 
the  transmitter  also  hops  L/2/10  times,  fin  our  case  there 
are  JO  bits  in  one  hop.)  The  receiver  still  have  to  catch 
up.  The  time  needed  is  L/2/10  bits.  Therefore,  the  total 
time  =  (L/2)+(L/2/J0)=fL/2) *1 .1  bits. 


Table  ]_ 


Variance  of  Acquisition  Time  (Experiment  Result) 
Hop  rate  =  30  hops/S,  2-tap  Model 


Efa/n  (db) 


Variance 


0.1 

0.418418 

0.02 

0.418418 

0.01 

0.303081 

0.005 

0.275799 

0.002 

0.435575 

0.1 

0.455856 

0.02 

0.337374 

0.01 

0.412955 

0.005 

0.921933 

0.002 

2.070953 

0.1 

0.524299 

0.02 

0.778812 

0.01 

1.200690 

0.5 

6.588859 

0.3 

5.536296 

0.1 

1.629105 

0.02 

6.980358 

Normalized 
Standard  Deviation 


0.62966 

0.62966 

0.55344 

0.55099 

0.56009 


0.61307 

0.50216 

0.55545 

0.62320 

0.51520 


0.65931 

0.69792 

0.63962 


0.93680 

0.97243 

0.82622 

0.95249 


0.5 

6.665808 

0.3 

2.210138 

0.1 

0.840163  . 

0.5 

9.171375 

0.3 

4.471619 

0.1 

7.753166 

0.5 

16.20247 

0.3 

13.74216 

0.1 

18.79460 

36.24931 

48.56720 

146.13722 


0.89386 

0.76763 

0.56669 


0.89559 

0.83897 

0.93249 


0.77684 

0.94278 

0.84126 


0.79713 

1.15495 

0.97152 


Table  2 


Variance  of  Acquisition  Time  (Experiment  Result) 
Hop  rate  =  100  hops/S,  4-tap  Model 


Eb/n  (db) 


Variance 


Normalized 
Standard  Deviation 


0.056443 


0.68907 


Tabl ' 


(Con  t ’ d  . ) 


iy  n  (db)  I 

PD  ; 

.  u.ce 

Normal i zed 

Land  a  rd  Deviation 

0.02 

0.028083 

0.55752 

30 

0.01 

0 . 039172 

0.60551 

0.005 

0.039606 

0.51557 

0.1 

0.022601 

0.56278 

25 

0.02 

0.033616 

\ 

0.52921 

.  0.01 

0.069021 

0.59486 

0.1 

0.045693 

0.62705 

‘  20 

0.02 

0.106365 

0.64079 

;  :  0.5 

0.308951 

0.86610 

15 

.  •  0.3 

0.289192 

1.00024 

:  0.1 

0.146960 

0.80117 

0.5 

. 

0.358209 

0.99068 

12 

.  0.3 

0.144358 

0.77637 

O.i 

0.201937 

0.69824 

0.5 

■  0.811281 

0.92320 

9 

0.3 

;  0.485079 

0.94808 

0.1 

0.388513 

0.55792 

0.5 

1.514624 

0.87441- 

6 

;  0.3 

1.176961 

1.00146 

0.1 

.  2.291038 

0.84144 

3 

0.5 

3.302562 

0.90081 

0.3 

2.015820 
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Fig.  1  The  receiver  in  FH  SS 
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Fig.  3  Hopping  frequency  slot 
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